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Flow Equations without Mean Field Ambiguity 
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We compare different methods used for non-perturbative calculations in strongly interacting 
fermionic systems. Mean field theory often shows a basic ambiguity related to the possibility to 
perform Fierz transformations. The results may then depend strongly on an unphysical parameter 
which reflects the choice of the mean field, thus limiting the reliability. This ambiguity is absent 
for Schwinger-Dyson equations or fermionic renormalization group equations. Also renormalization 
group equations in a partially bosonized setting can overcome the Fierz ambiguity if the truncation 
is chosen appropriately. This is reassuring since the partially bosonized renormalization group ap- 
proach constitutes a very promising basis for the explicit treatment of condensates and spontaneous 
symmetry breaking even for situations where the bosonic correlation length is large. 

PACS numbers: ll.10.-z, ll.10.Hi, ll.10.St 



I. INTRODUCTION 

Mean field theory is a widely used method in statis- 
tical physics and quantum field theory, in particular for 
ground states characterized by condensates and sponta- 
neous symmetry breaking. For example, mean field so- 
lutions of the Nambu-Jona-Lasinio (NJL) model Q or 
extensions of it are one of the main theoretical tools in 
nuclear physics. The recent discussion of color supercon- 
ductivity at high but realistic baryon density is mainly 
based on this method 0, 0, 0, IS H • One out of many 
examples from statistical physics is a mean field descrip- 
tion y| of antiferromagnetic and superconducting con- 
densates in the Hubbard model 

II Hi EH Quite gen- 
erally, mean field theory (MFT) seems to be well suited 
for systems with multifermion interactions and bosonic 
condensates. 

It is well known that MFT has a basic ambiguity which 
is connected with the possibility to perform Fierz trans- 
formations (FT) for the underlying local multifermion 
interaction. For example, in the Hubbard model this 
ambiguity has a sizable influence on the results p|. The 
origin of this ambiguity becomes apparent already in the 
simplest NJL-type model (for only one fermion species) 
with a chirally invariant pointlike four fermion interac- 
tion: 



S F = / d A x 



1 



V^# + ^A CT [(^) 2 -(^) 21 



(1) 



Due to the Fierz identity 

[WW) 2 - (vvVV) 2 ] +2 [(^) 2 - W^V) 2 ] = o (2) 

only two of the quartic couplings arc independent and we 
write 

Act = X a + 27 Ay, Ay = (l-7)Ay, Xa = 7Ay- (3) 

The parameter 7 is redundant since it multiplies just the 
vanishing expression J5J. No physical quantity can de- 
pend on 7 in a full computation of the functional integral 
for the partition function and expectation values of field 
operators. The model is completely characterized by the 
two "physical" couplings A CT and Ay. 

We will see in sect. [H] that the MFT results can 
strongly depend on 7, limiting their quantitative re- 
liability. For example, the value of the critical cou- 
pling A™* for the onset of a nonvanishing condensate 

cr ~ (qf) ^ 1 ~ 7 ^ tp) depends on 7 for fixed Ay ^ (cf. the 

values in the first row of Tabs. ITT HT|l . This "Fierz ambi- 
guity" is a strong effect unless |Ay| is much smaller than 
I Act I . Similarly, for given Ay, X a the mean field value of 
the condensate a in the phase with spontaneous symme- 
try breaking depends on 7. 

The origin of the mean field ambiguity is easy to un- 
derstand: once the pointlike interaction has been "dis- 
tributed" on the channels (S - P), (V), (A) (with re- 
spective couplings A CT , Ay, Xa) the interactions in the 
channels (V) and (A) do not influence 1 the computation 
in the channel (S — P) anymore - at least as long as the 



In this paper we concentrate on this model which is rcg 
ularized by a sharp momentum cutoff q 2 < A 2 . 
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1 We refer here to MFT as used in most computations and im- 
plemented on a more formal level by the Hubbard-Stratonovich 
transformation or partial bosonization. Sometimes the wording 
"mean field" is also used for a Schwingcr-Dyson approach for 
which no ambiguity is present. We discuss this in sect. IVIIII We 
note that the MF-ambiguity gets enhanced once we include, in 
addition, di-fermion channels ~ ipip and i/>-0. 
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mean fields ~ ip'y^t/j and ^^^ip do not yet get an ex- 
pectation value. This distribution depends, however, on 
the parameter 7. It may sometimes be possible to de- 
velop from other considerations an educated guess what 
should be the distribution of the interaction on the var- 
ious channels, thus limiting the range of "acceptable" 7. 
Still, the spread of the results over the acceptable range 
of 7 should be considered as a lower bound for the sys- 
tematic uncertainty of a MFT computation. Depending 
on the other uncertainties the relative importance of this 
"Fierz- ambiguity" may be more or less important. For 
the example of Tab. [I]the Fierz ambiguity of MFT seems 
to be of the same size as the spread in the results between 
the different methods beyond MFT. On the other hand, 
for the values used in Tab. [H]thc ambiguity for the range 
< 7 < 1 is so large that no quantitative statement is 
possible for MFT unless 7 can be restricted to a much 
smaller range. 

A correct "guess" of 7 is often difficult 2 . One may 
therefore prefer to reduce the dependence on the unphys- 
ical parameter 7 by computing in a more elaborate ap- 
proximation. This situation is very similar to perturba- 
tive computations in quantum field theory: the results in 
a given order depend on the renormalization scheme and 
the renormalization scale [i. Typically, the scheme de- 
pendence (/x-dependence) gets reduced in higher orders. 
Since the physical results cannot depend on /1 the residual 
/i-dependence is often used as a guess for the remaining 
error - in fact it constitutes a lower bound for the sys- 
tematic uncertainty in a given order of the perturbative 
expansion. 

In this paper we want to investigate methods where 
MFT appears as some type of first step in a more system- 
atic expansion. We find indeed that for these methods 
the dependence on 7 is reduced as compared to MFT. For 
some approximations a residual 7-dependence remains 
which should again be considered as a lower bound on the 
systematic uncertainties at this level 3 . We will discuss 
methods based on the exact renormalization group equa- 
tion for the effective average action 01 or 011 Schwinger- 
Dyson equations |l3l Il4j . To improve the approximation 
we include additional diagrams similar to those needed 
at order ^ in the ^-expansion [Tl liH liH HI Il9j . 

As a guidance, we first study in sect. IIVI perturbation 
theory in the one-loop approximation. While the results 
are independent of 7, the validity of the perturbative cal- 
culation is limited to small coupling. In particular, the 
interesting phenomenon of spontaneous symmetry break- 



ing (SSB) cannot be seen. 

This shortcoming is improved substantially by the 
use of a non-perturbative flow equation for the scale- 
dependence of the effective average action 0|. In 
this approach an infrared cutoff with scale k is intro- 
duced such that only fluctuations with momenta q 2 > k 2 
contribute effectively. The effective couplings therefore 
depend on the scale and their fc-dependence is governed 
by a renormalization group (RG) equation. For example, 
in a simple truncation takes the same form as the ac- 
tion P) ; but the couplings become now scale dependent, 
i.e. \ a — ► A CT f. etc.. In the limit k — the infrared cutoff 
is absent and all fluctuations are included. For k = the 
effective average action Tk=o is the generating functional 
of the lPI-Grcens functions of the full theory. 

The scale dependence of Tk is described by an exact 
functional differential equation |l2j |. Non-perturbative 
approximations to this exact equation involve a trunca- 
tion of the general form of Tk- The non-perturbative 
flow equations for the quartic couplings A CTl fc, Xv,k and 
\a,u ar e solved in sect. |V] Now the onset of SSB is sig- 
naled by a divergence of the quartic couplings A. The 
fermionic flow equation yields results for the critical cou- 
plings which are independent of 7. Unfortunately, in con- 
trast to MFT, one cannot easily extend the investigation 
to the phase with SSB and compute, for example, the 
value of the order parameter or the fermion mass gap. 
This would require the inclusion of multifermion interac- 
tions in the flow (e.g. eight fermion interactions), leading 
to high algebraic complexity. 

Partial bosonization seems to be the ideal remedy to 
this difficulty [H il IH HI |1 IH. In sect. |mj the 
model JTJ is rewritten as an equivalent Yukawa type 
model with scalars <j>, vectors V M and axial vectors 
representing the corresponding fermion bilincars. Spon- 
taneous symmetry breaking can now be dealt with by 
computing the effective potential for <f> and looking for a 
minimum at <f> ^ 0. For example, a term ~ <p stands 
for an eight fermion interaction. Unfortunately, par- 
tial bosonization brings back the "Fierz ambiguity" of 
MFT. In fact, an approximation which only includes the 
fermionic fluctuations and omits the bosonic fluctuations 
is precisely equivalent to MFT. Nevertheless, MFT can 
now be considered as a starting point of a more system- 
atic procedure which also includes the bosonic fluctua- 
tions. Here again, a test for the validity of any approxi- 
mation to the bosonic fluctuations should see that phys- 
ical quantities become independent of 7. 

In sects. IVII IVIII we study the flow equations for the 
Yukawa model with action 



2 If the effective four fermion interaction is known beyond the 
pointlike limit the distribution on the various channels is much 
more restricted. It should be approximated by a sum of terms 
where Xy only depends on the total momentum of the bilinear 
■07 M i/) and similar for the other channels. 

3 Other uncertainties not related to the Fierz ambiguity may still 
be larger. This is obviously the case for methods without a Fierz 
ambiguity. 
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this model is equivalent to the NJL-type model QJ. 
We choose a simple truncation for I\. where only the 
(fc-dependent) couplings appearing in Eq. 10} are re- 
tained. In sect. IVII we first neglect the fact that new 
quartic fermion interactions are generated by the flow. 
Including the bosonic fluctuations leads to a running of 
the Yukawa couplings and reduces substantially the 
dependence of the results on 7 as compared to MFT. Still, 
the flow equations are not one-loop exact in this trunca- 
tion, as also reflected by the residual 7-dependence (cf. 
the third row in Tabs. [I][nJ). The reason for the incom- 
pleteness at the one-loop level is the omission of quartic 
fermion cou pling s which are generated by the flow for 
k < A. In [2y| a systematic method has been devel- 
oped regarding how these interactions can be eliminated 
by fc-dependent field redefinitions, leading to a modifica- 
tion of the flow equation. We use this method in sect. 
IVIII Indeed, the modified flow equations in the trunca- 
tion Q are now one-loop exact. They are equivalent to 
the fermionic flow equation of sect. At this point we 
have reached a satisfactory starting point for an exten- 
sion of the flow equation beyond the four fermion trun- 
cation of sect. Ivl and beyond the MFT- type truncation 
of sect. IVIII Further extensions of the truncation in the 
bosonic sector are now straightforward and will be briefly 
discussed in sect. IIXI They are, however, not the main 
focus of this paper. 

We investigate in sect. IVIIII bosonization from the 
viewpoint of the Schwinger-Dyson equations 0, IT3 |. 
Those, too, can be applied to both the purely fermionic 
and the partially bosonized setting. The formulation us- 
ing the fermionic model £Q| is naturally unambiguous 
even in the simplest approximation. In the bosonic model 
Q we investigate two very simple approximations, the 
simplest of which corresponds to MFT. The other one 
takes into account mass corrections due to diagrams with 
internal boson lines. In a very simple approximation for 
the bosonic propagator and vertices it corresponds to the 
gap equation derived for the fermionic model. The re- 
sults are therefore independent of 7. Finally, we present 
in sect. |X]an overview over the merits and shortcomings 
of the different methods. 



II. CRITICAL COUPLINGS FROM MEAN 
FIELD THEORY 

A mean field calculation treats the fermionic fluctua- 
tions in a homogenous background of fermion bilinears 

£ = $ & = -® {}¥-) V>>, % = (ihrf) 

and = (V>7 AI 7 5, 0)- It seems straightforward to replace 
in the four fermion interaction in Eq. Q one factor by 



(■0'0) 2 — ($j 5 ip) 2 -> 2H>(1 +7 5 )?A- 2^-0(1 -7 5 )V', 

$7 M 7fy) 2 -» 2A^ 7 5 ip. (6) 
The partition function becomes then a functional of 0, 

Z$, V,A] = J Vi>V^ exp (-S$, V, 4>, V, i]) , (7) 

where S is given by Eq. with the replacements 10 • 
Self-consistency for the expectation values of the fermion 
bilinears requires 

4> = ^(l-T^^A-^lnZ, (8) 

and similar for the other bilinear A^. Chiral symmetry 
breaking by a nonzero 4> requires that the " field equation" 
© has a nontrivial solution. We note that Z[<j>, V, A] cor- 
responds to a one-loop expression for the fermionic fluc- 



(F) 



tuations in a bosonic background. With T\ 
the field equation is equivalent to an extremum of 



AnZ 



(F) 



r( F ) = J d 4 x h\,p$ + -\ v v lt v» + -x A A^j+r\ 

(9) 

A discussion of spontaneous symmetry breaking in MFT 
amounts therefore to a calculation of the minima fr( p ). 

We note that this calculation can be done equivalcntly 
in the Yukawa theory Q, JSJ. The mapping of the 
bosonic fields reads cf) = {h a / y* 2 ) <j>, = (/iy/^y)V^, 

A^ = {Ha I Ma) -^m- Keeping the bosonic fields fixed and 
performing the remaining Gaussian fermionic functional 
integral yield precisely Eq. © . Mean field theory there- 
fore corresponds precisely to an evaluation of the effective 
action in the partially bosonized Yukawa model in a limit 
where the bosonic fluctuations are neglected. 

We want to compute here the critical couplings (more 
precisely, the critical line in the plane of couplings A CT , 
Ay) for which a nonzero expectation value </> ^ indi- 
cates the onset of spontaneous symmetry breaking. For 
this purpose we calculate the mass term ~ </>*</> in r^ F ' 
and look when it turns negative. This defines the crit- 
ical couplings. We assume here a situation where the 
expectation values of other bosonic fields such as V M or 
A M vanish in the relevant range of couplings. It is then 
sufficient to evaluate r^ F ' for V M = A^ = 0. 

In a diagrammatic language gaussian integration over 
the fermionic variables corresponds to evaluating the dia- 
gram of Fig. 13 We define our model with a fixed ultravi- 
olet momentum cutoff q 2 < A 2 , such that the MFT-result 
becomes (1*4 = l/(327r 2 ), x = q 2 ): 
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FIG. 1: Bosonic mass correction due to fermion fluctuations. 
Fermionic lines are solid with an arrow, bosonic or "mean field 

lines" ~ ( i^.) i>) are dashed. 



r A2 

/ da;a;lii(a; + /i^». (10) 
Jo 



From this one finds the mean field effective action 

= rfVrf 5 (11) 

= (/4 - Avihlk 2 ) cf>*ct> + const + O ((</>» 2 ) , 

where we have expanded in powers of <fi. The mass term 
turns negative if 



2/4 
KIN 2 



< 8v 4 , 



(12) 



As it should be this result only depends on the ratio 
hl/nl = 2A CT . 

We now want to determine the critical line in the plane 
of invariant couplings \ a , Ay from the condition 1)12(1. i.e. 



A crit = 



8w 4 A 2 ' 



(13) 



Using the relation we infer a linear dependence of 
Acrit Qn ^ whenever Ay ^ 



1 



8w 4 A 2 



- 2 7 Ay. 



(14) 



(For numerical values see Tabs. HI and HT|) . This depen- 
dence is a major shortcoming of MFT. We will refer to it 
as the "Fierz ambiguity". The Fierz ambiguity does not 
only affect the critical couplings but also influences the 
values of masses, effective couplings etc.. 

The origin of the Fierz ambiguity can be traced back 
to the treatment of fluctuations. A FT of the type 
changes the effective mean field. In a symbolic language 
a FT maps (i) a i/> a )(i>bipb) — * {^a^i^bipa) where the 
brackets denote contraction over spinor indices and ma- 
trices ~ 7 M or r~j 7 5 are omitted. A mean field 4> a ipa-, 
appears after the FT as i/j a ijjb- From the viewpoint of the 
fluctuations one integrates out different fluctuating fields 
before and after the FT. It is therefore no surprise that 
all results depend on 7. 



III. PARTIAL BOSONIZATION 

The MFT calculation introduces "mean fields" com- 
posed of fermion - antifcrmion (or fermion - fermion) bi- 
lincars. This is motivated by the fact that in many physi- 
cal systems the fermions are not the only relevant degrees 
of freedom at low energies. Bosonic bound states become 
important and may condense. Examples are Cooper pairs 
in superconductivity or mesons in QCD. For a detailed 
description of the interplay between fermionic and com- 
posite bosonic fluctuations it seems appropriate to treat 
both on equal footing by introducing explicit fields for 
the relevant composite bosons. This will also shed more 
light on the status of MFT. 

Partial bosonization [2(| El HU is achieved 4 by 
introducing unit factors in the functional integral for the 
partition function 



Z = 



VipV<ipexp(-S[ip]) 



(15) 



= J VipVipV<j)VV> 1 VA l *Af ct> Af v Af A exp(-S[iP}) 



with 



A/0 = exp 



JVv =exp 



+ 2^(1 + 7 5 )V; 
- 2^(1 -7 5 )V> 



(16) 



x V, 



Ha = exp 



/V 

4 (a* - %^7 5 ^ 

V Va 

The action in the " bosonic language" is composed of the 
original fermionic action and the arguments of the expo- 
nentials in A/0, Ay and Ma- Using the relation JSJ) one 
finds that the quartic fermion interaction is cancelled. It 
is now replaced by mass terms for the bosons and Yukawa 
couplings between bosons and fermions as given by the 
expression Q. 

Since the action, Eq. is quadratic in the fermionic 
fields the functional integral over the fermionic degrees 
of freedom is Gaussian and can be done in one step. As 
we have seen in the previous section this precisely leads 
to the MFT results. More precisely, we understand now 
that for different choices of 7 the MFT treatment leaves 
out different bosonic fluctuations. In this context we note 



4 The inclusion of fermionic and bosonic source terms as well as 
infrared cutoffs is straightforward and omitted here for simplicity. 
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that the bosonic potential in the bosonized effective ac- 
tion must be bounded from below. This restricts the 
possible couplings to A CT , Ay , A^ > 0. In the invariant 
variables this restriction translates to X a , Ay > and for 
7 it implies < 7 < 1. 



IV. PERTURBATION THEORY 

In order to cure the unpleasant dependence of the MFT 
result on 7 wc will include part of the bosonic fluctua- 
tions in sects. IVII and lYlII Some guidance for the level of 
approximations needed can be gained from perturbation 
theory in the fcrmionic language. Since the four fermion 
vertex is uniquely characterized by X a and Ay the pertur- 
bative result must be independent of 7 at any given loop 
order. The lowest order corrections to the four fermion 
couplings are obtained by expanding the one-loop expres- 
sion 5 



Ar (i-ioo P ) = _ STr 
2 



lnfS< 2 > 



-Tr 



(17) 



up to order (ipip) 2 . The corresponding graphs are shown 
in the second panel of Fig. From 

Ar (i-ioo P ) = W4A a| ( 18 ) 

[4A 2 - 4\ a (\ A - 2Ay)] [0W>) 2 ~ tyftf 
+ [-2A CT Ay + A(X A - X v )Xv] [(h^) 2 
+ [-A 2 + 2X a X A + 3X 2 V - 2X A X V - A|] [(^W*) 



we can read off the corrections AA^ , AAy and AA^ to the 
coupling constants. In order to establish that our result 
is independent of 7 we use the freedom of FT to bring 
our results into a standard form, such that AAa = -2—. 
Inserting next the invariant variables leads to: 



AA CT = 4v 4 A 2 (A 2 +4A a Ay + 3A y ), 
AAy = 2u 4 A 2 (A (T + Ay) 2 . 



(19) 



In contrast to MFT the result does not depend on 7. 

The perturbative result, Eq. flTfy . always leads to fi- 
nite corrections to the coupling constants. Remember- 
ing that in the fermionic language the onset of SSB is 
marked by a divergence of the coupling constants, it be- 
comes clear that we will never get SSB in perturbation 



5 The supertrace STr pro vide s an appropriate minus sign in the 
fermionic sector (see e.g. 1271 1. Moreover, in the full S^ 2 ^ matrix 

wc have a term from the 5 2 /S'ip5tp derivative (S^,) and a term 
from the 8 2 /Stf>8ip, accounting for a factor of 2 in the language 
with the normal trace. The trace includes momentum integration 
and summation over internal indices. 



theory. No critical couplings can be calculated. This is a 
severe shortcoming of perturbation theory which cannot 
be overcome by calculating higher loop orders. Only an 
infinite number of loops can give SSB. In the next section 
we establish how a renormalization group treatment can 
overcome this difficulty without encountering the Fierz 
ambiguity of MFT. A calculation of the critical coupling 
becomes feasible. Nevertheless, even this RG-treatment 
has its limitations once the couplings diverge. In par- 
ticular, it does not allow us to penetrate the phase with 
SSB. In sects. I VII and TVIII this shortcoming will be cured 
by a RG-treatment in the partially bosonized language. 
In particular, we will see in sect. IVIII which diagrams are 
needed in order to maintain the independence of results 
on 7 in analogy to perturbation theory. 



V. RENORMALIZATION GROUP FOR 
FERMIONIC INTERACTIONS 

The renormalization group equations for the effective 
average action [l2| are obtained by adding a fc-dependent 
infrared cutoff term AS k to the action. The effective 
avera ge a ction Tk results 6 from a Legendre transform of 
\nZk Due to the infrared cutoff, receives only 

contributions from fluctuations with q 2 > k 2 . We take 
ASfc quadratic in the fermion fields, 



AS* 



(2tt)' 



■$(q)Rk(q)ip{q) 



(20) 



with Rk vanishing fast for q 2 k 2 . For a cutoff function 
which diverges for k — > 00 and vanishes for k — > the 
functional Tk interpolates between the "classical action" 
Too = Sf and the full effective action To = T which in- 
cludes all quantum fluctuations. The fc-dependence of Tfc 
is governed by the following exact equation (t = ln(fc)): 

dtTkbP] = ^STr{d t \n(T^^}+R k )}, (21) 



d t = d t R k 



d 
dR k 



(2) 

with T k the second functional derivative of T k with re- 
spect to the fermion fields. Despite its suggestive one- 
loop form this is a functional differential equation which 
cannot be solved exactly. 

A first approximation neglects the fc-dependence of T k 
on the right hand side (RHS). The solution is the pertur- 
bative result, Eq. (|17|l . As we have seen in the previous 
sect. II VI this approximation does not lead to SSB. For a 
better approximation we restrict T k to the terms specified 
in Eq. Q but take all couplings explicitly fc-dependent. 



6 We do not differentiate here the notation between the arguments 
of Tfc and the fluctuating fields. In terms of the arguments of T k 
the cutoff is substracted in the definition of Y k Il2j . 
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In the action Q we have only local interactions. Ex- 
pressed in momentum space the four fermion interactions 
have no momentum dependence. This is often referred to 
as the local potential approximation (LPA) |2£| |2^, |3(j ■ 

(2) 

Decomposing the fluctuation matrix T k according to 

if ' + R k = V + T (22) 

into a field independent part V (inverse propagator) and 
a field dependent part T we can expand the RHS of Eq. 
(|2*T|) as follows: 

d t r k = \sTr{d t ^\ -\sTr{d t {^A } (23) 

This amounts to an expansion in powers of fields and we 
can compare the coefficients of the four fermion terms 
with the couplings specified by Eq. JQ). We obtain a set 
of ordinary differential equations for the couplings: 

d t ~Xv,k = -Av i l[ F)A (s)k 2 C\^ k + X v ,k) 2 , (24) 

in agreement with j3l"l | where the same model has been 

studied. The threshold functions l^' 4 are defined in |3^ . 
For our actual calculation we use a linear cutoff 7 [13 and 
adapt the threshold functions to our setting with fixed 
momentum cutoff q 2 < A 2 in App. [X] The dependence 
on s = k 2 1 A 2 becomes relevant only for k > A whereas 
for k < A one has constants l[ F ^ 4 = 1/2. 

The fermionic flow equations 8 (|24() do not depend on 
7. In a diagrammatic language we again have evalu- 
ated the diagrams of Fig. 0] (second panel) but now with 
fc-dependent vertices. In the RG-formulation we only go 
a tiny step Afc, and reinsert the resulting couplings (one- 
loop diagrams) before we go the next step. This leads to 
a resummation of loops. Since Eq. I|24|) is now nonlin- 
ear (quadratic terms on the RHS) the couplings can and 
do diverge for a finite k if the initial couplings are large 
enough. Therefore we observe the onset of SSB and find 
a critical coupling. Since Eq. I|24|) is invariant this criti- 
cal coupling does not depend on 7! Values for the critical 
coupling obtained by numerically solving Eq. i|24|) can be 
found in Tabs. [Qand[n] 



7 The threshold functions depend on the precise choice of the cut- 
off. For the very simple truncation used in this paper this de- 
pendence can actually be absorbed by a suitable rescaling of k, 
cf. App. EI 

8 As discussed above, the perturbative result, Eq. 1191 . can be 
recovered from Eq. 1241 if we neglect the fc-dependence of the 
couplings on the RHS and perform the t-intcgration. 



The next step in improving this calculation in the 
fermionic language would be to take the momentum de- 
pendence of the couplings into account (e.g. j^il) or to 
include higher orders of the fermionic fields into the trun- 
cation. This seems quite complicated and at first sight 
we have no physical guess what is relevant. The renor- 
malization group treatment of the bosonic formulation in 
sect. IIIII seems much more promising in this respect. 

VI. BOSONIC FLOW 

The flow equations in the bosonic language are ob- 
tained in complete analogy with the fermionic formula- 
tion. In this paper we restrict the discussion to a " point- 
like" truncation as given by Eq. with fc-dependent 
couplings. We will see that in this approximation we 
reproduce the result of the last section if we take care 
of the fact that new fermionic interactions are generated 
by the flow and have to be absorbed by an appropri- 
ate fc-dependent redefinition of the bosonic fields. The 
crucial advantage of the bosonic formulation is that it 
can easily be extended. For example, the bosonic bound 
states become dynamical fields if we allow for appropriate 
kinetic terms in the truncation, i.e. 

Ar t ,„=y <f Jz^* a<v + ^fv r ,.v»* + ^a^a<™ 

with 

V ia , = d v y v -d v V Vki = d M A v - d v A„. (26) 

Also spontaneous symmetry breaking can be explic- 
itly studied if we replace /i 2 0*0 by an effective poten- 
tial C/((^*</>) which may have a minimum for 0^0. 
This ap proach has been followed in previous studies 
[35L l36t l37l l3S[ . In the present paper we will briefly 
pursue only the extension (|25|1 since our main focus is 
the issue of the Fierz ambiguity. 

It is instructive to neglect in a first step all bosonic 
fluctuations by putting all bosonic entries in the propa- 
gator matrix V~ x equal to zero. This removes all dia- 
grams with internal bosonic lines. Among other things 
this neglects the vertex correction, Fig. [21 an d therefore 
the running of the Yukawa couplings. Indeed, Fig. ^is 
the only diagram contributing and we recover MFT. One 
obtains the flow equations 

d t ^ k = 8h 2 aM , A k 2 l[ F ^(s), 

d t ^ k = 8h 2 VM v i k 2 l[ F) '\s), 

dtA, k = 8hX kVi k 2 l{ F)A ( s ), 

& t K, k = 0, d t h v ,k=0, d t h A , k =Q. (27) 

As long as we do not consider the wave function renor- 
malization l|25|l for the bosons, the flow can be completely 
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Using the dimcnsionlcss e's wc now find: 



FIG. 2: Vertex correction diagram in the bosonized model. 
Solid lines are fermions, dashed lines are bosons. There exist 
several diagrams of this type since we have different species 
of bosons. 



described in terms of the dimcnsionlcss combinations 



ev,k 



/4,fc 

h 2 k 2 ' 
n V,k K 



£A,k 



f*A,k 



(28) 



Due to the constant Yukawa couplings we can integrate 
Eq. (|27|) . We find critical couplings: 



hlA 2 



h 2 v A 2 



(29) 

These are, of course, the results of MFT, Eq. We 
note that in Eq. I|27|l the equations for the different 
species of bosons are completely decoupled. The mass 
terms do not turn negative at the same scale for the 
different species. Indeed it is possible that the mass of 
one boson species turns negative while the others do not. 
Such a behavior is expected for the full theory, whereas 
for the fermionic RG of sect. |V]all couplings diverge si- 
multaneously due to their mutual coupling. However, no 
real conclusion can be taken from Eq. (|29|l because of 
the strong dependence on 7. 

Now, let us take also the bosonic fluctuations into ac- 
count. This includes the vertex correction, Fig. [21 and 
the flow of the Yukawa couplings docs not vanish any- 
more. In the pointlike approximation (Z$ = Za = Zy = 
0) one obtains 



dth% k = 



d t h 2 Vk 



- 32v 4 l\ h {s)k 2 h 



d t h 2 



t"-A,k 




(30) 



Qt^a.k = 

d t £v,k = 
dtf-A.k = 



-2e CT , fc + 8 
-2ev !k + 8 
-2i A ,k + 8 



1 + 4 

2 + 

2 - 



ey,fc eU.fc 
t-v.k , ey,fc 

£A,k 
£A,k 

ev,k 



2e<T,fc 

^A,k 



2r 



v4 F) '\s). 
(31) 



The onset of spontaneous symmetry breaking is indicated 
by a vanishing of e for at least one species of bosons. 
Large e means that the corresponding bosonic species 
becomes very massive and therefore effectively drops out 
of the flow. 

For initial couplings larger than the critical values (see 
Tabs. [I] and [HJ both and ey,k reach zero for finite 
t. Due to the coupling between the different channels 
they reach zero at the same t. At this point eA.k reaches 
infinity and drops out of the flow. This is quite different 
from the flow without the bosonic fluctuations where the 
flow equations for the different species were decoupled. 
The breakdown of all equations at one point resembles 9 
now the case of the fermionic model discussed in sect. 
IV! The 7-dcpendence of the critical couplings is reduced 
considerably, as compared to MFT. This shows that the 
inclusion of the bosonic fluctuations is crucial for any 
quantitatively reliable result. Nevertheless, the difference 
between the bosonic and the fermionic flow remains of the 
order of 10%. 



VII. ADAPTED BOSONIC FLOW 

In our truncation the bosonic propagators are approx- 
imated by constants /i^ 2 . The exchange of bosons there- 
fore produces effective pointlike four fermion interactions. 
One therefore would suspect that this approximation 
should contain the same information as the fermionic 
formulation with pointlike four fermion interactions. An 
inspection of the results in Tabs. [n] shows, however, 
that this is not the case for the formulation of the pre- 
ceding section. In particular, in contrast to the fermionic 
language the results of the bosonic flow equations still 
depend on the unphysical parameter 7. 

In fact, even for small couplings A the bosonic flow 
equations of sect. do not reproduce the perturba- 
tive result. The reason is that at the one-loop level new 
quartic fermion interactions are generated by the box di- 
agrams shown in Fig. |3J An easy inspection shows that 
they contribute in the same order A 2 as the diagrams in 
Figs. ^ and [21 Even if we start from vanishing quartic 
couplings after partial bosonization, such couplings are 



This is an artefact of the pointlike approximation. 
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(a) 



(b) 



FIG. 3: Box diagrams for the bosonized model. Again, 
solid lines are fermions, dashed lines bosons and vertices are 
marked with a dot. The diagrams generate new four fermion 
interactions even for the model @ without direct four fermion 
interactions. 



generated by the flow. The diagrams in Fig. [3] yield 



dtK.k = Px„ = -%v A l[ F) '\s)e a 2 K +4fc-*7(fc), 

Mo-.fe ^A,k 



h 2 h 2 

2 'V,k ,l A,k 



d t X A k 



2 2 
f^V.k l L AM 



(32) 



t4,k 



h 4 h 4 
12-^ - 12- A 



f*v,k 



+ 2fc~ 2 7(fc). 



Here j(k) is an in principle arbitrary function of scale de- 
termining the choice of FT for the generated four fermion 
interactions. We will make a special choice of this func- 
tion (similar to the one made in sects. IIVI and lV|) namely 
we require 



7 



1-7 



(33) 



with e given in Eq. 128fl . The resulting equation 
dt(iv,k/^A,k) = fixes 7(fc). An improved choice of 7(fc) 
can be obtained once the momentum dependence of ver- 
tices is considered more carefully |26| . 

An inclusion of the couplings X k into the truncation of 
the effective average action does not seem very attractive. 
Despite the partial bosonization we would still have to 
deal with the multi-fcrmion interactions and the bosonic 
formulation would be of even higher algebraic complexity 
than the fermionic formulation. A way out of this has 
been proposed in |26|. There, it has been shown that 
it is possible to reabsorb all four fermion interactions 
generated during the flow by a redefinition of the bosonic 
fields. In the following brief description of this method 
we use a very symbolic notation. 

Introducing an explicit fc-dependence for the definition 
of the bosonic fields in terms of fermion bilinears, the flow 
equation Eq. 1)2 111 is modified: 



d t r k = d t r k 



ocpk 



(34) 



Here dtT k \,p k = dtT k | is the standard flow of the effective 
average action at fixed fields. Shifting <f> by 

d t 4>k = (i'lf) d t uj k (35) 

we find 

cV=d tM 2 |, (36) 
d t h = d t h | +fi 2 d t u>k, d t X = d t X | -hd t oj k 
and we can choose co k to establish: 

d t \ = 0. (37) 

Instead of including running four fermion couplings ex- 
plicitly we therefore only have to use adapted flow equa- 
tions for the couplings contained in Eq. |0J. 

Let us now apply this method explicitly to our model. 
Shifting 

d t (j> = -if) ( - X ^ if)d t uj at k, (38) 



d t <j>* = i/> 



1 — 7' 
_ 2~ 

1+7 5 



h„_ k d t uJa.k-, (39) 



we have 

OtXa,k = 9(A CTl fe| — n a> kOt^(7,k, 

dtXy,k = d t Xv,k\ — 2hy,kdt0Jv,k, 
d t X A M = d t X A ,k\ ~ 2h A , k d t ujAM- 

Requiring dtX = for all A's we can determine the func- 
tions u>: 

P\ a a . Av Q . Px 

h a ,k 



d t UJa„ 



d t uJv,i 



9 t WA,k 



(40) 



2hv,k 2h A ,k 

with the /?- functions given in Eq. (|32|l . This yields the 
adapted flow equations for the Yukawa couplings 

dtK.k = d t h a ,k\ + /4,fcdtkV,fc, (41) 
d t hv,k = d t h v ,k\ + (J-v^dtujv^k, 
d t h Atk = d t h Atk \ + [i\ k d t io A ,k- 
Combining Eqs. (2D), (E2J), (E3- E3, 6U determines 

m 

j(k) = 2 Vi l[ F) -\s) (42) 



1 



~Vk e V,k€A,k 



(4g CT;fc - i Aik ) 2 

£A,k(ev,k + £A,k) 



Having fixed the ratio between ey.fc and e A ^ k we need 
only two equations to describe the flow. We will use the 
ones for l a , k and l v .k = (1 - l)^v,k 

d t l a , k = ~2e a , k + 4[(l + 7 ) -4(-2 + 7 + 2 7 2 )^ 
+ 4(3- 77 + 4 7 3 )|^]4 F) ' 4 ( S )«4, 



d t £v,i 



c V,k 

+ 4[^-(2 7 -i)] 2 /r ),4 (^4. 



2l„ k 



(43) 
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These equations are completely equivalent to the 
fermionic flow Eq. I|24|) . In order to see this we remem- 
ber that the simple truncation of the form J3J is at most 
quadratic in the bosonic fields. We can therefore eas- 
ily solve the bosonic field equations as a functional of 
the fermion fields. Reinserting the solution into the ef- 
fective average action we obtain the form Q with the 
fc-dependent quartic couplings 



W,k 



k 2 ev,k ' 



(44) 



Inserting this into Eq. (|43|) we find Ea . (|24|) . establishing 
both the exact equivalence to the fermionic model and 
the 7-indcpcndcncc of physical quantities. 

On this level of truncation the equivalence between the 
fermionic and the adapted bosonic flow can also be seen 
on a diagrammatic level. As long as we do not have 
a kinetic term for the bosons the internal bosonic lines 
shrink to points. On the one-loop level we find an exact 
correspondence between the diagrams for the bosonized 
and the purely fermionic model summarized in Fig. 0] 
This demonstrates again that one-loop accuracy cannot 
be obtained without adaptation of the flow. 



VIII. THE SCHWINGER-DYSON APPROACH 

Finally, we compare in this section the results of MFT 
and the RG-treatment with the Schwinger-Dyson (SD) 
approach [Tsl IT^j . On the exact level the RG and SD 
approaches are equivalent in the sense that the propa- 
gator and higher N-point functions calculated using the 
flow equation J^TJ) are also solutions of the SD-equations 
[39L l40| . Nevertheless, once truncations are used the re- 
sults will, in general, differ. The RG-cquations resum 
diagrams beyond the leading order SD-cquation. 

We start with the purely fermionic model For 
this model the Schwinger-Dyson equation, approximated 
to lowest order, is depicted in Fig. 



5(a) 



It is a closed 

equation since only the bare four fermion vertex appears. 
(Only higher order terms involve the full four fermion 
vertex.) We write the full fermionic propagator Gf as 



G? 1 b)=Gp 1 (p) + S F (p) 



(45) 



with the free propagator Gfo and self-energy Ep- Using 
this one obtains a gap equation for the self-energy which 
can be solved self-consistcntly. To simplify the discussion 
we make an ansatz for the self-energy: 



M F 7 5 



(46) 



where the effective fermion mass M-p obeys the gap equa- 
tion 



f A2 

AIf = 8U4 [A<j + Ay] / dx : 
Jo 



■ M * (47) 
x + M 2 ' 



The onset for nontrivial solutions determines the critical 
couplings: 



[\ a + \v] ( 



1 



8v 4 A 2 ' 



(48) 



This result is shown in Tabs. HI ITU and does not depend 
on 7, as expected for a fermionic calculation. We observe 
that the MFT-result for the A" u coincides with the SD- 
approach for a particular choice 7 = 1/2. However, in 
general MFT is not equivalent to the lowest order SD- 
equation. This can be seen by computing also the criti- 
cal coupling for the onset of SSB in the vector channel. 
The MFT and SD results do not coincide for the choice 
7=1/2. 

Next, we turn to the SD-equations for the bosonized 
model Q which are depicted in Fig. 5(b) We will make 



here two further approximations by replacing in the last 
graph of Fig. 



5TbJ 



the full fcrmion-fcrmion-boson ver- 
tex by the classical Yukawa coupling and the full bosonic 
propagator by ^ig 2 . We remain with two coupled equa- 
tions. 

In a first step we approximate these equations even 
further by neglecting the last diagram in Fig. 5(b) al- 
together. Then no fermionic propagator appears on the 
right hand side of the equation for the fermionic propa- 
gator which only receives a mass correction for (</>) 7^ 0. 
Without loss of generality we take </> real such that 
Mp = h a (j) and 



G F -\q) = -{l + h a( h 5 - 



(49) 



Inserting this into the equation for the expectation value 
6 we find 



4^4 

14 



dxx- 







+ hi <b 2 ■ 



(50) 



For the onset of nontrivial solutions we now find the crit- 
ical value 



Wo 



1 



crit 



8 Vi A 2 



\X a + 2 7 A y 



crit 



(51) 



which is the (ambiguous) result from MFT given in Eqs. 
(UJJ and (HJ. This is not 

surprising since this exactly is 
MFT from the viewpoint of Schwinger-Dyson equations. 
Indeed, Eq. (|50|) is precisely the field equation which 
follows by differentiation of the MFT effective action © 
with respect to <fi, 



r< F ) = nl<f> 2 - Avi / dxx\n{x + h 2 J 2 ). 



(52) 



In a next step we improve our approximation to include 
the full set of diagrams shown in Fig. 5(b) Using the 



same ansatz as before the self-energy £p now has two 
contributions, 



S F = A/ F 7 5 = h a (f>Y + Am F 7 



(53) 
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FIG. 4: Summary of all diagrams encountered in the previous sections. There is a one to one correspondence between the 
diagrams of the bosonized model (first row) and the purely fermionic model (second row). Solid lines with an arrow denote 
fermionic lines. The letters in the diagrams are given to visualize the ways in which the fermionic operators are contracted, 
e.g. the first diagram in the second row results from a term [(•0 a , a )('i/) c i/) c )][(i/) c i/) c )('i/)i,i/)f,)]. Shrinking bosonic lines (dashed) to 
points maps the diagrams in the first row to the second row. In the approximations of sect. IVII only the first or the first two 
diagrams are taken into account. 




FIG. 5: Diagrammatic representation of the lowest order Schwinger-Dyson equations for the fermionic model (a) (Eq. (0) 
and the partially bosonized model (b) (Eq. |QJ). The shaded circles depict the full propagator, the circle with the cross is the 
expectation value of the bosonic field and the empty circle is the full Yukawa vertex. 



The first one is the contribution due to the expectation 
value of the bosonic field whereas Amp is the contribu- 
tion from the last diagram in Fig. |5(b) given by an 
integral which depends on Mp. Both in the equation for 
(4>) and in the equation for the fermionic propagator only 
Mp appears on the RHS. Inserting (4>) in the graph, Fig. 
|5(b)| one finds a gap equation which determines Mp: 



Mp = &v 4 



3. 

2/4 



My 



Ma 



dxx- 



Mp 



(54) 



Once more, this can be expressed in terms of the invariant 
couplings and again we arrive at Eq. 147II . 



Looking more closely at the two contributions to Mp 
we find that alone neither the contribution ~ <fi (which 
amounts to MFT as we have discussed above) nor the 
"fermionic contribution" Amp is invariant under FT's. 
Only the combination Mp, which is the fermion mass 
and therefore a physical quantity, is invariant. Indeed, 
changing the FT amounts to a redefinition of the bosonic 
fields. In a very similar fashion to what we have done in 
sect. IVIII it allows us to choose bosonic fields such that 
Amp = 0. Taking 7=1/2 gives us such a choice of the 
bosonic fields. This explains why MFT gives the same 
result as the purely fermionic calculation in this special 
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case. 



IX. BEYOND THE POINTLIKE 
APPROXIMATION 

In order to compare the size of the Fierz ambiguity 
with typical errors from the truncation we extend in this 
section our investigation beyond the pointlikc approxima- 
tion. We add to the general form of the effective action 
the bosonic propagator terms (|25(1 with av = a a = 1. In 
order to simplify the numerical computations we choose 
in this section the ERGE-regularization scheme (cf. App. 
|B)| rather than a sharp ultraviolet momentum cutoff. As 
discussed in App. [51 the definition of the "classical cou- 
plings" A(j, Xy differs from the sharp cutoff rcgulariza- 
tion. 

We have extended the truncation in consecutive steps 
in order to make the effect of various truncations di- 
rectly visible. Our results are shown in Tab. UnJ By 
comparing the first four rows of Tab. IIIII with Tab. [HJ 
we observe the effect of the different ultraviolet regu- 
larization. With the ERGE-dcfinition of the "classical 
couplings" one finds a substantially larger value of the 
critical coupling A™* than for the sharp cutoff rcgular- 
ization. This different size is purely a matter of defini- 
tions and not related to any approximation error. Let 
us now extend the truncation in the partially bosonized 
formulation step by step. As discussed above, the lowest 
order is MFT. As a typical size for the Fierz ambiguity 
we quote A\ = A™* (7 = 0.75) - A= rit (7 = 0.25) = -20. 
The next step (1) includes the bosonic fluctuations in the 
pointlike approximation, without the adaption discussed 
in sect. IV11I At this level the residual Fierz ambiguity is 
substantially reduced, Aa = 1.71. Step (2) includes the 
adaption (box diagrams) in the pointlike approximation. 
Due to the equivalence with the pointlike approximation 
in the purely fermionic langauge there is no Fierz ambi- 
guity at this level. 

As a first step beyond the pointlike approximation we 
include the running 10 of the wave function renormaliza- 
tions (WFR) Z$, Zy, Za- (Details of the flow equations 
will be published elsewhere.) In step (3) this is done with- 
out the adaption by "rebosonization" . Omitting com- 
pletely the rebosonization ((l)+(3)) yields A\ = 1.57 
whereas combining Eqs. (1), (2) and (3) results in 
A\ = —0.31. We finally use rebosonization also for the 
momentum dependence in the four fermion and Yukawa 
couplings generated by the flow (box diagrams and vertex 
corrections). The results depend somewhat on the de- 
tailed method (to be published elsewhere) and one finds 
A a = 0.22. It is impressive how the truncation reduces 
the Fierz ambiguity from a value Aa = —20 for MFT to 



For simplicity we neglect the anomalous dimensions in the 
threshold functions. 



a value |A A | < 1! 

On the other hand, a comparison of the last rows in 
Tab. 11111 yields a typical value for the remaining trunca- 
tion error, A" 1 ' = 58 ± 5. We emphasize that the quoted 
error should not be taken as an (more sophisticated) error 
estimate in a strict sense. A true error estimate is noto- 
riously difficult in a situation without small parameters. 
One possibility would be further extensions of the trunca- 
tion. One could also investigate the influence of different 
rebosonization procedures or the choice of different cutoff 
functions Rk(q)- Nevertheless, it seems convincing that 
we have reached a truncation uncertainty that is far less 
than the Fierz ambiguity in MFT (i.e. \A\\ = 20). On 
the other hand, on this higher level of the truncation the 
residual Fierz ambiguity is already much less (|Aa| < 1) 
than the truncation error. This is precisely what one 
would like to achieve for a successful non-perturbativc 
approximation scheme in a partially bosonized setting. 

In passing, we note that Schwinger-Dyson equations 
do a pretty good job in our context. 



X. SUMMARY AND CONCLUSIONS 

We compare different approximation methods for a 
strongly interacting fermion system - the NJL model in 
our case. For this purpose we have computed the criti- 
cal couplings for the onset of chiral symmetry breaking 
using mean field theory (MFT), fermionic and bosonic 
renormalization group methods (RG) and the Schwinger- 
Dyson equation (SD). This permits a direct compari- 
son between the various methods. We believe that the 
general characteristics found here remain valid for other 
strongly interacting systems as well. For a sharp momen- 
tum cutoff the results for A" 1 * are summarized in Tabs. 
[I][n]for two fixed values of Ay. Corresponding results for 
the ERGE rcgularization can be found in Tab. IIIII Since 
the most characteristic features and problems of the dif- 
ferent methods are most clearly seen when the couplings 
A CT and Xy are of similar size we concentrate the discus- 
sion on Tab. [HJ 

All methods discussed here (except perturbation the- 
ory in sect. HVfl correspond to non-perturbative resum- 
mations of perturbative diagrams. Both MFT and the 
lowest order SD sum only over fermionic fluctuations in 
presence of a bosonic background. They include, in prin- 
ciple, the same type of diagrams, Fig. The MFT- 
result depends strongly on the choice of the background 
field. This "Fierz ambiguity" is expressed by the depen- 
dence on the unphysical parameter 7 in the tables. No 
such ambiguity appears in the SD approach which there- 
fore seems, at least at first sight, more reliable. We note 
that for a particular choice of 7 the MFT and the SD ap- 
proaches give identical results - in our case 7 = 1/2. This 
has led to widespread belief that MFT and SD are equiv- 
alent if the basis for the Fierz ordering is appropriately 
chosen. However, this is not the case, as can be seen by 
calculating also the critical coupling where spontaneous 
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Approximation 


Sect. 


7 = 


0.25 


0.5 


0.75 


1 


MFT 


LLU 


39.48 


38.48 


37.48 


36.48 


35.48 


Ferm. RG 


El 


41.54 


41.54 


41.54 


41.54 


41.54 


Bos. RG 


|VT| 


36.83 


36.88 


36.95 


37.02 


37.12 


Adapted Bos. RG 


IVIU 


41.54 


41.54 


41.54 


41.54 


41.54 


SD 


Ivnii 


37.48 


37.48 


37.48 


37.48 


37.48 



TABLE I: Critical values A" 1 * for Xv = 2 and for various 
values of the unphysical parameter 7 (with A = 1). Progress- 
ing from MFT to the bosonic RG and adapted bosonic RG 
the dependence on 7 decreases as more and more diagrams 
are included. The Schwinger-Dyson result is independent of 
7 but contains no vertex corrections in contrast to the Re- 
calculations. 



Approximation 


Sect. 


7 = 


0.25 


0.5 


0.75 


1 


MFT 


LLU 


39.48 


29.48 


19.48 


9.48 


-0.52 


Ferm. RG 


El 


14.62 


14.62 


14.62 


14.62 


14.62 


Bos. RG 


|VT| 


15.44 


13.39 


13.45 


15.55 


19.46 


Adapted Bos. RG 


IVIU 


14.62 


14.62 


14.62 


14.62 


14.62 


SD 


Ivnii 


19.48 


19.48 


19.48 


19.48 


19.48 



TABLE II: The same" as in Tab. Qbut with X v = 20. 

"The negative sign for the critical coupling at 7 = 1 in the MFT 
calculation means that the system is in the broken phase for any 
positive value of A CT in this calculation. 



symmetry breaking sets in in the vector channel (in the 
absence of other order parameters). There is again a 
value 7 = — ( A(j + Ay)/(2Ay) where MFT and SD give 
identical results, but it differs from 7 = 1/2 as encoun- 
tered in the scalar channel 11 . We conclude that there is 
no possible choice of 7 where both critical couplings for 
SSB in the scalar and vector channels are identical in the 
MFT and SD approaches. The conceptual and practi- 
cal difference between the two approaches appears even 
more clearly if we consider a model with eight-fermion- 
couplings instead of a quartic coupling. Whereas in MFT 
the onset of SSB can be computed in one-loop order, 
only a three loop diagram contributes to the gap equa- 
tion for My in the SD approach. We also warn that the 
choice of 7 for which MFT and SD coincide in a given 
channel is not necessarily the optimal choice. Comparing 
the SD result (or the MFT result for 7 = 1/2), namely 
A (j = 19.48, with the perhaps more precise result from 
the renormalization group, X a = 14.62, we see that an 
"optimal choice" of the FT for the MFT approach may 
correspond to a value of 7 above 1/2. 

Partial bosonization is a very powerful tool for under- 
standing strongly interacting fermionic systems beyond 
the level of MFT or SD-equations. It allows us to treat 
the bosonic fluctuations in an explicit manner and pro- 
vides for a rather simple framework for the discussion of 



11 Actually, 7 is negative and therefore outside the range of strict 
validity of MFT. 



Approximation 


Chap. 


7 = 0.1 


0.25 


0.5 


0.75 


0.9 


MFT 


mi 


74.96 


68.96 


58.96 


48.96 


42.96 


SD 


|V111| 


58.96 


58.96 


58.96 


58.96 


58.96 


Bos. RG=(1) 


EH 


53.16 


52.93 


53.32 


54.64 


55.88 


Ad. Bos. RG=(l) + (2) 


El 


58.83 


58.83 


58.83 


58.83 


58.83 


(l) + (3) 


lixj 


53.90 


53.66 


54.00 


55.23 


56.37 


(l)-(3) 


lixl 


58.14 


58.04 


57.88 


57.73 


57.64 


(1H4) 


iLXl 


61.60 


61.69 


61.82 


61.91 


61.94 



TABLE III: Critical coupling X" lt for Xv = 20 for a UV reg- 
ularization by the ERGE scheme We show different approxi- 
mation steps: (1) the pointlike contibutions to the mass and 
the Yukawa coupling (Figs. (2) the pointlike contribu- 

tions from the box diagrams (Fig. |3J, (3) the contribution 
to the WFR from the purely bosonic diagram (Fig. [TJ and 
(4) the contribution to the WFR from the momentum depen- 
dence of the diagrams |21 and |3] (those contribute again via an 
appropriate adaption of the flow). 



SSB. Most importantly, it permits the direct exploration 
of the ordered phase which is, in practice, almost inac- 
cessible for the fermionic RG. In order to permit a simple 
comparison with the fermionic RG we have used a very 
crude approximation for the purely bosonic sector by re- 
taining only a mass term and neglecting bosonic interac- 
tions as well as the momentum dependence of the bosonic 
propagator. In this approximation the effect of the boson 
exchange between fermions does not go beyond pointlike 
fermionic interactions. Taking into account only the run- 
ning of the Yukawa couplings (Fig. |2J in the bosonic RG 
of sect. IVII we observe already a very substantial de- 
crease of the Fierz ambiguity as compared to MFT. The 
dependence on 7 is greatly reduced and the numerical 
value of the critical coupling comes already close to the 
result of the fermionic RG. These features can be com- 
pared to the inclusion of higher loop effects in pertur- 
bation theory in particle physics: they often reduce the 
dependence of the results on unphysical parameters, such 
as the choice of the renormalization scale. 

As compared to perturbation theory, the box diagrams 
© are still missing in the discussion of sect. IVII This 
shortcoming is cured by the adapted bosonic renormal- 
ization group discussed in sect. IVIU Here the relation 
between the bosonic composite fields and the fermion bi- 
linears becomes scale dependent. This formulation is well 
adapted to the basic idea of renormalization where only 
effective degrees of freedom at a certain scale k and their 
effective couplings should matter for physics associated 
with momenta q 2 < k 2 . The system should loose all 
memory of the detailed microscopic physics. In partic- 
ular, the choice of an optimal bosonic field for the long 
distance physics should not involve the parameters of the 
microscopic theory, but rather the renormalized parame- 
ters at the scale k. In this formulation it has also become 
apparent that the distinction between "fundamental de- 
grees of freedom" and "bound states" becomes a mat- 
ter of scale j2(|. The adapted bosonic RG reproduces 
in our crude approximation the results of the fermionic 
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RG. Wc argue that for precision estimates in the partially 
bosonized approach the " adaptation" of the definition of 
the composite field seems mandatory. 

It seems plausible to us that within the local interac- 
tion approximation considered in this paper the most re- 
liable results are obtained from the fermionic or adapted 
bosonic RG of sects. IVland lVIll First, these methods 
sum over a larger class of diagrams. The diagrams in- 
cluded by the other approaches are all contained in the 
ones taken into account by the fermionic or adapted 
bosonic RG. Second, the renormalization procedure ac- 
counts properly for the fact that the relevant physics de- 
pends on scale. In the low momentum region relevant for 
spontaneous symmetry breaking all physics should be de- 
scribable in terms of effective low energy couplings. This 
is not realized by the MFT or SD approaches where the 
microscopic or "bare" couplings appear explicitly. We 
point out, however, that beyond the local interaction ap- 
proximation the relative merits of the various methods 
discussed here depend on the physical situation of the 
investigated model. 

The aim of the investigation of the bosonic RG in this 
paper is, of course, not simply a reproduction of the re- 
sults of the fermionic RG. By a careful comparison be- 
tween different approaches we rather want to open the 
door for future more elaborate techniques for the study of 
strongly interacting fermion systems. With the present 
results the adapted non-perturbative flow equations in 
the partially bosonized approach offer an ideal start- 
ing point for an investigation of spontaneous symmetry 
breaking. Without too much effort we can now include 
the momentum dependence of the bosonic propagator 
which goes beyond the approximation of local fermionic 
interactions. This is crucial for the understanding of crit- 
ical phenomena for which the renormalized boson mass 
vanishes. We have made a first step in this direction 
by including the wave function renormalization in sect. 
IIXI Furthermore, bosonic interactions can now be in- 
cluded in the form of an effective potential for the scalar 
field. This is mandatory for the RG-investigation of the 
ordered phase where the minimum of the potential oc- 
curs for <fi ^ 0. Both effects have already been treated 
for the NJL-model with three colors and two flavors for 
the fermions (quarks) [25t 135^ . In the view of the re- 
sults of the present paper, the quantitative accuracy and 
conceptual setting of such investigations could further be 
improved by the " adaptation" of the effective composite 
variables. Including these effects one by one will also 
permit a more detailed appreciation of the uncertainties 
remaining in the present truncation of the flow. 

In summary, some of the methods proposed to deal 
with strongly interacting fermionic systems have an ad- 
ditional source of systematic uncertainty. The Fierz am- 
biguity is related to the choice of the bosonic or mean 
field, parametrized by an unphysical parameter 7. For 
those methods the spread of the results within an ac- 
ceptable range of 7 should be considered as a lower bound 
for this additional systematic uncertainty. Different val- 



ues of the parameter 7 correspond to an identical initial 
fermionic action. Therefore, vanishing or at least small- 
ness of the Fierz ambiguity should be required for the 
self-consistency of an approximation. We find that, de- 
pending on the model and parameters, MFT can have a 
very substantial ambiguity which should then be reduced 
by systematic improvements. 

On the other hand, the Fierz ambiguity is, of course, 
not the only source of error - several methods as SD or the 
fermionic RG have no such ambiguity by construction. 
Without a systematic error analysis which is highly diffi- 
cult for non-perturbative systems and beyond the scope 
of this paper there is no simple overall criterion in order 
to judge which method is most suitable - the answer may 
actually depend on the detailed model and problem. The 
non-perturbative flow equations based on an exact renor- 
malization group equation offer interesting prospects for 
an understanding of problems where a large correlation 
length is important. It is reassuring that the analysis 
of the potential Fierz ambiguity in a partially bosonized 
setting shows that this approach can be considered as a 
promising extension beyond MFT. 
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Note added 

A renormalization group approach close to the spirit of 
the SD-equations has been proposed recently in the 
context of the " bosonic effective action" . There the long 
range bosonic fluctuations can be incorporated without 
invoking partial bosonization. The Fierz ambiguity is 
therefore absent in this approach. 



APPENDIX A: THRESHOLD FUNCTIONS FOR 
FINITE UV-CUTOFF 

In this appendix we compute the threshold functions 
as defined in for the linear cutoff given in [3i| in the 
presence of a finite UV-cutoff A. The inverse (massless) 
average propagator Pb for bosons and the corresponding 
squared quantity Pp for fermions are given by 

Pb = q 2 + Z-lR k (q)=q 2 (l + r B (q 2 )), (Al) 
P F = q 2 (l + r F (q 2 )) 2 , 

where tb and rp reflect the presence of the IR-cutoff. 
The dimensionlcss functions tb and rp only depend on 
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y = q jk . For the linear cutoff [33j they read 



r B (q,k) = (--1)6(1- y), 



(A2) 



r F (q,k) = [ — -l)e(l-y). 



In presence of an ultraviolet cutoff A and in the absence 
of mass terms the threshold functions can only depend 
on the ratio s = k 2 /A 2 . With 



q 2 d[Z,j, tk r B ] d 



2 P F d[Z^ ik r F \ d 



Z^ k dt dP B Z 4 , M l + r F 
we find for bosons {x = q 2 ), 

j(B)d , \ 



dt 



dP F 
(A3) 



(A4) 



= h 
2 

_ 2 

~~ d 



dxx%- 1 dt\n(P B {x) +u)k 2 ) 



d + 2 



1 



1+UJ 



-9(1 



(2 + djl-s- 1 ))^ 
2(d+2) 



1 



1+UJ 



e( fl -i) 



and for fcrmions, 



,(F)d 



(A5) 



2 



dxx?- 1 d t hi(P F {x) + ujk 2 ) 



1 



1 + UJ 



0(1 -a) 



(d + 1 - rfs- 1 )^ 



1 



1 



1 + w 



6(s-l) 



Higher threshold functions can be obtained simply by 
differentiating with respect to ui: 



1 



n + 8 n o da; 



^(w,t?.«). 



(A6) 



For a finite value of the UV-cutoff A the threshold func- 
tions are explicitly s- and therefore fc-dependent. Taking 
A — > oo wc have s = for any value of fc. This renders the 
threshold functions k- independent. In the present work 
we neglect the anomalous dimensions 77^, 77^, and effec- 
tively only consider a fermionic cutoff since Z^.k = 0. For 
our purpose the fermions are massless and we abbreviate 
for lu = 0, 



#(0,0, a) =#(*). 
This yields explicitly 



(A7) 



l[ F) ' i (s) = l[Q(l-s)+s- 2 Q(s-l)]. (A8) 



To obtain the pcrturbative result from the fermionic RG- 
equation we used 



dtk 2 l[ F) 



(s) = / dkkl[ F)A {s) = 



A2 
2 ' 



(A9) 



As long as we keep the sharp momentum cutoff at q = A 
this integral is universal, i.e. it does not depend on the 
precise choice of the IR-cutoff. Indeed the universality 
is necessary to reproduce perturbation theory for every 
choice of the IR-cutoff. 

We have also used other cutoff functions Rk different 
from the linear cutoff. Within the local interaction ap- 
proximation we have found that the value of the critical 
coupling comes out independent of the choice of Rk ■ The 
basic reason is that a multiplicative change of l[ F ^' 4 due 
to the use of another threshold function can be compen- 
sated by a rescaling of k (cf. Eq. I|24|ll. The rescaling 
is simply multiplicative for s < 1. with a suitable gener- 
alization for s > 1. Critical values of the flow which are 
defined for k — > 00 are not affected by the rescaling. Let 
us demonstrate this for Ay . Writing Eq. (|2"4*f) in the scale 
variable k we have 



d k X v ,k = -±v A l[ F)A {s)k{\ a , k + \ v ,kf 



Rescaling to 



fc(fc) 



dk kl 



(f),4 



s) 



wc find 



d k X V,k 



(A10) 



(All) 



(A12) 



Due to the universality of Eq. (|A9|I the domain for k, 
[0,oo], is now mapped to [0, ^-], giving the domain for 
k independent of the IR-cutoff. Having obtained identi- 
cal differential equations for every choice of the IR-cutoff 
without any rescaling of A establishes the above claim for 
the critical couplings. 

Note however, that this would not hold if we would 
start the integration of the flow equation at k = A. In this 
case the domain [0, A] for k is mapped into an interval for 
k that depends on the threshold function and therefore on 
Rk ■ Actually, the Rk dependence in this case is not very 
surprising because different IR-cutoffs then correspond 
to different UV regularizations (see App. 0. Since our 
model is naively non-renormalizable results can depend 
on the choice of UV regularization. 



APPENDIX B: UV REGULARIZATION 
ERGE SCHEME 



THE 



In the first eight sections we have implemented the UV 
regularization by a sharp cutoff in all integrals over mo- 
mentum space. This is often not the most practical reg- 
ularization. As an alternative, the modes with q 2 > A 2 
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FIG. 6: In a UV regularized theory not all modes contribute 
completely. This plot schematically depicts "how much" each 
mode contributes. The thick line is for the sharp momen- 
tum cutoff. All modes with q 2 < A 2 are included completely. 
Other UV regularizations (dashed, thin dashed and thin solid 
line) typically not only include a small fraction of the high 
momentum modes, q 2 > A 2 , but in addition leave out a small 
fraction of the low momentum modes. 



are not completely left out, only suppressed, as depicted 
in Fig. El 

This can be realized by a modification of the propaga- 
tor in the action yielding an "U V regularized classical ac- 
tion" . Different UV regularizations usually correspond to 
different "classical actions" . Therefore, it is no surprise, 
that different UV regularizations give different results. 
In particular, this is true for the critical coupling in the 
NJL model, Eq. (JTJ, as one can see by comparing Tabs. 
H] [H] calculated using a sharp UV cutoff at A with Tab. 
IIIII which employs a UV regularization by the ERGE 
scheme (below) with the linear cutoff Eq. 1A2(I . One 
expects that dimcnsionlcss low energy quantities such as 
the ratio of particle masses and order parameters in the 
phase with spontaneous symmetry breaking, show much 
less dependence on the regularization pBj . 

In App. |J|we have evaluated the threshold functions 
for a theory which is UV regularized by a sharp cutoff 



in momentum space. The threshold functions depend 
on the ratio s = p- in a rather complicated way. Cutoff 
scale independent threshold functions would be desirable, 
to simplify numerical calculations. 

This is implemented by the ERGE regularization 
[42l . l43l | where A is taken to infinity in the definition of the 
threshold functions (i.e. s = 0). On the other hand, the 
classical action is now specified indirectly by the "initial 
value" of the effective average action T\ at some ultra- 
violet scale A. Depending on the choice of the cutoff 
function Rk(q) the fluctuations with momenta q 2 > A 2 
have not yet been integrated out completely, as depicted 
in Fig. As a result, the relation between "classical 
couplings" A CT , Xy and physical observables depends on 
the choice of Rk and differs from the regularization with 
a sharp cutoff. We emphasize that this is a difference in 
the definition of the classical couplings and should not 
be confounded with approximations errors. 

Although it is usually not the simplest method, we can 
invoke UV regularization by the ERGE scheme also in 
the context of perturbation theory or SDE. This follows 
along the lines indicated in Sect. [V]for perturbation the- 
ory. Typically any expression can be written in terms of 
inverse propagators P, internal momenta q we integrate 
over, and external momenta p we do not integrate over, 



F(P,q,p). 



(Bl) 



A specific ERGE scheme is specified by the choice of 
the IR regulator Rk- Replacing the inverse propaga- 
tor P by the IR regularized inverse propagator P + Rk 
we can calculate the contribution from each scale fc, 
k~ 1 dtF(P + Rk,q,p). Integrating over all scales from 
ko = A to k = we obtain the UV regularized expres- 
sion, 



dkk^dt 



k =A 



/ F(P 



Rk,q,p) 



(B2) 
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